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Abstract 

The flow into an open return wind tunnel 
inlet was simulated using Euler equations An 
explicit predictor-corrector method was employed 
to solve the system. The calculation is time- 
accurate and was performed to achieve a steady- 
state solution. The predictions are in reasonable 
agreement with the experimental data Wall 
pressures are accurately predicted except in a 
region of recirculating flow. Flow-Held sur- 
veys agree qualitatively with laser velocimctcr 
measurements The method can be used in 
the design process for open-return wind tun- 
nels 

Introduction 

The design of low-speed wind tunuels has 
received a great deal of attention recently because 
of the need for improved and larger wind tun- 
nels to advance the technology of low-speed flight 
In the past, the design process has been more of 
an experimental procedure, which requires exten- 
sive model and prototype testing Computational 
fluid dynamics, however, is now sufficiently ad- 
vanced so it can be used to assist such a design 
process and thereby reduce the amount of testing 
required The purpose of the present investiga- 
tion is to numerically analyze the aerodynamics of 
the inlet of the 80x120 open return wind tunnel 
of the National Full-Scale Aerodynamics Complex 
(NFAC) as part of an effort to improve the test- 
section flow quality. 

Since the aerodynamics of wind tunnels is com- 
plex, it is divided into a study of the components; 
for example, into the individual studies of inlet, 
contraction nozzle, diffuser, etc. The inlet flow 
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here constitutes the flow in the inlet section, con- 
traction nozzle, and test-section of a How-through 
facility. The purpose of the inlet is to take air 
from the atmosphere, condition it for the wind 
tunnel circuit, and accelerate it through the bell- 
mouth to the test-section with minimal turbulence 
and minimal large scale disturbances The inlet 
flow is therefore strongly dependent on the na- 
ture of the atmospheric winds and turbulence, and 
on the inlet geometry The downstream effect 
of the fan and diffuser section on the inlet How 
will usually be negligible The test-section flow 
quality is evaluated in terms of the flow steadiness, 
uniformity, angularity, turbulence level, boundary 
layer growth, etc . The present numerical study 
was undertaken to examine the ability of an Euler 
method to accurately predict some of these fea- 
tures in the flow Held of the inlet part of the cir- 
cuit 

Problem Deflmtion 

The inlet flow is simulated using the unsteady 
compressible Euler equations in three dimensions. 
The solution is obtained by marching in tune 
to a steady state. The equations are sohed 
using MacCormack’s explicit predictor-corrector 
scheme 1 in Cartesian coordinates with a finite 
volume formulation. The Euler code 2 is rnodiDcd 
to admit a zonal calculation without having to 
read in from or out to a disk This is made 
possible by defining a single index for all the 
grid points in the computational domain, thereby 
reducing the storage requirement The computa- 
tions were carried out on the Cray-XMP com- 
puter at NASA Ames Research Center The 
code was used to simulate the flow induced in 
the wind tunnel by a fan downstream of the 
test section (sink-type flow). The predictions 
are compared with' the experimental data, the 
pressures and the velocities are compared at 
different locations, and the velocity vectors are 
plotted to compare the predicted flow pattern 
with that obtained by laser velocimeter measure- 
1 ments. 



Formulation 

The three-dimensional, unsteady Euler equa- 
tions in conservation form and in cartesian coor- 
dinates are given by. 


d t Q + d x E + d y G + d,H = 0, (1) 
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Q is the solution vector, E, G and H are the fluxes 
in the x, y and z directions respectively, and the 
pressure is given by 


the unit vector normal to the area element 
dS and S is the surface bounding the volume 
V. 

Predictor-Corrector Scheme 

The finite-volume, explicit-integration scheme 
used in the present study is described in Ref. 
2. The basic numerical algorithm is due to 
MacCormack. 1 The flux integral in equation (3) 
is approximated by summing the dot products of 
the area vectors and the appropriate flux vectors 
over each of the three orthogonal projections of 
the volume element V Integration in time is 
done using a two-step predictor-corrector sequence 
Although the formulation is time-accurate, only 
steady-state solutions are sought The predic- 
tor advances the solution to time level n -J- 1 
as 



and the corrector updates this solution as 


P = (l~ l)(e- 
where q 2 = u 2 -f- v 2 -f- w 2 

The physical variables above are normalized with 
respect to the free-stream density and free-stream 
speed of sound c Writing equation (1) as 

d t Q-\-V F = 0 (2) 

where F — ( E,G,H ), the Euler equations can 
be expressed in integral form by integrating 
them over a small stationary volume element V 
as 

Sv d ‘ Q + lv V - F = 0 

Using the mean value and divergence theorems, the 
integral form reduces to 

d t Q = -~fn-FdS (3) 

V J S 

where Q is redefined to be the mean value 
over the volume element V and is evaluated 
at some point interior to V (assumed to be 
its centroid for all practical purposes), H is 
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where F is the flux vector evaluated at the 
"centroid” of the volume element V and F+ is 
the flux vector evaluated at the centroid of the 
volume element next to V m the increasing coor- 
dinate direction given by t. Similarly, F~~ cor- 
responds to the volume element next to V in the 
decreasing coordinate direction; S is the area vec- 
tor of any side of the hexahedral volume element in 
question The subscript p refers to the values at the 
predicted level. The time step criterion for At is 
given by the Courant-Friednch-Lewy stability con- 
dition. The maximum allowable time step is given 

by 


I?-5| + C|S| 

where 5 is the velocity vector 

Grid Generation 

The three-dimensional grid for the inlet, Fig 
1, is generated algebraically The grid (shown in 
Fig 2) represents the ground plane, and the top 
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and side computational boundaries upstream of the 
inlet, the symmetry plane midway between the side 
walls of the tunnel, the inlet cowls in the horizontal 
and the vertical directions, the nozzle and the test- 
section. The nozzle has a 5 1 contraction ratio 
Two planar grids are also shown in Figs 3 and 
4. The z-plane represents a horizontal cut, and 
the y-plane represents a v ertical cut, Fig 3(b) and 
Fig 4(b) show enlarged views near the cowl The 
grid is stretched so as to resolve the flow gradients 
in more detail The three-dimensional grid is 
generated plane by plane using an interpolation 
scheme. 

Initial Conditions 

The computations of the inlet sink-type flow 
are carried out corresponding to quiescent condi- 
tions at the upstream boundary Initially, stag- 
nation conditions exist everywhere. The pres- 
sure drop which generates the flow in the inlet 
and yields a given test-section velocity is intro- 
duced at the downstream boundary gradually in 
time The flow can be imagined to exit at the 
downstream boundary into ambient conditions at 
somewhat lower than the upstream stagnation pres- 
sure and is given by p e = p t — J pq e 2 in the 
steady incompressible limit, where the subscripts 
e and s refer to the exit and stagnation condi- 
tions 

Boundary Conditions 

The boundary conditions at the wall are imple- 
mented using a dummy volume concept The 
condition of surface tangency at solid and sym- 
metry boundaries is enforced using the vector equa- 
tion 


{pq)o = (p?)» - 2 
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where the subscripts o and i refer to the outer 
(dummy) and inner volume elements at the boun- 
dary, ^ is the velocity lector and 5 is the area 
vector at this boundary across which the density, 
pressure and energy are reflected 2 This condi- 
tion reflects the mass flux vector across the boun- 
dary. 

The implementation of the solid boundary 
condition at the convex corners (e g., B and C 
in Figs 3 and 4 respectively) is done in such a 
way that the same image (dummy) volume ele- 
ment is used to satisfy the boundary conditions 
(corresponding to surface tangency) at the two in- 
terior nodes adjacent to the corner This is done se- 
quentially in the three coordinate directions as the 
fluxes are accumulated in the integration process, 


thereby removing the corner discontinuity. The 
integration process is similar to that used in Ref 
3, except that here all the fluxes are accumu- 
lated prior to updating Around such regions as 
the external cowl, where the curvature is rela- 
tively large, a simple reflection boundary condi- 
tion on pressure 2 is no longer valid. Therefore, 
the pressure at the image node next to the wall 
is found from the conservation of momentum nor- 
mal to the cowl surface, which reduces to a balance 
between the centrifugal force and the pressure 
gradient in that direction. 4 This condition is given 
by 


9nP = 
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which is written in difference form as 


Po = P. + An 


pV 2 


for concave surfaces and as 


Po — Pi - An 


P V 2 


for convex surfaces, where n is the direction nor- 
mal to the surface, An is the distance between 
interior and dummy volume element centroids in 
the normal direction, V is the tangential velocity 
at the curved surface, and r is the radius of 
curvature of the surface. For the downstream 
boundary, a nonreflecting boundary condition for 
the pressure 5 was used while the other four 
independent variables were simply extrapolated 
The boundary condition on pressure is given 
by 


d t p — pcd t u + a{p — p e ) = 0 

where a is a parameter which determines the rate 
at which the transients die out before steady state 
is attained. 5 

The split-flux characteristic boundary condi- 
tions 6 were used at the inflow boundary. Stagnation 
pressure, stagnation enthalpy, and the flow direc- 
tion are prescribed Pressure is related to velocity 
along characteristics with a locally isentropic flow 
assumption. The inflow boundary is placed 
sufficiently far upstream so that the velocity 
gradients are small there This permits a stretched 
grid at the inflow boundary. 
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Results and Discussion 

Computations were carried out at the test 
section Mach number, M = 0 15 The results 
from the computations are shown in Figs 5- 
11 The converged results correspond to a 
relative convergence criterion on pressure such 
that 


(p m+l_ p m) 


< 0 0001 


slowly decelerates as it enters the test section The 
tendency of the flow to be directed toward the 
center of the plane has also been observed in flow- 
visualization studies. The crossflow plane passing 
through the inlet cowl shows the flow following the 
contours of the cowl as expected A comparison 
between the crossflow plane velocity vectors at x — 
— 27ft from laser velocimeter measurements and 
the computations is shown in Fig 13 The predic- 
tions are in qualitative agreement with the experi- 
ment 

Conclusions 


over the entire flow field, where m is the time- 
step index. The velocity vector plots in the 
x-y and x-z planes as shown in Fig 5 and 
Fig. 6 indicate the correct qualitative behavior 
of the flow The separation bubble (shown in 
Fig. 5(c)) immediately downstream of the cowl, 
as observed experimentally, is not predicted in 
the computations since the effects of viscosity are 
not simulated Figure 7 shows a variation of 
the axial velocity normalized by the test-section 
velocity with the spanwise coordinate, y, at an 
x station situated 27 ft upstream of the inlet 
and close to the 2 = 0 plane The agreement 
with the experiments in this case is quite good 
Figure 8 shows a variation of the normalized axial 
velocity with the vertical coordinate z at x = 
—27 ft, near the symmetry plane between the 
two side walls of the tunnel. The predictions are 
again in good agreement with the experimental 
data. 

Tunnel floor pressure distribution and resul- 
tant velocity along the center line given by C p ( C p 
is defined with respect to the velocity and pres- 
sure at the test section) and V (at the test section), 
respectively, are plotted versus the downstream 
distance, x, from the cowl to the test section in 
Fig 9. The agreement between the present cal- 
culations and the calculations made using a panel 
method is good except in the region of recir- 
culating flow immediately downstream of the cowl. 
Pressure variation with the downstream distance 
(C p distribution) on the side wall centerline near 
the 2 = 0 plane is shown in Fig 10 The com- 
parison between the predictions and the experiment 
is good except near the separation bubble where 
the pressures are overpredicted relative to the ex- 
periment. Figure 11 shows the crossflow plane grid 
at three streamwise stations, x = 36, 130 and 
260 ft respectively Velocity vector plots at the 
corresponding streamwise stations in Fig 12 show 
the flow tending toward the center of the crossflow 
plane Figure 12 also shows that the crossflow ac- 
celerates in passing through the nozzle and then 


A three-dimensional, unsteady, inviscid flow 
code has been used to predict the pressures and 
the velocity profiles inside the inlet of the NFAC. 
The predictions are in reasonable agreement with 
the experiment. Although the viscous effects are 
not accounted for, many features of the inlet flow 
are captured satisfactorily The general flow direc- 
tion is predicted in accord with laser-doppler an- 
nemometer studies The pressure variation on 
the solid walls as well as the flow velocities are 
close to the experimental values However, vis- 
cous effects such as separated flow regions, secon- 
dary flows, and turbulence can be simulated only 
through a viscous flow computation However, 
using the present approach, various design criteria 
such as pressure losses and the desired test sec- 
tion flow uniformity can be estimated, given the 
geometry and the fan characteristics of an open- 
return circuit The results of the computations in- 
dicate that the Euler calculations can be used to 
simulate the flow in open-return wind tunnel in- 
lets 
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a) View showing the inlet model components: inlet cowl, 
nozzle, and the test section. 



b) View looking at the front into the tunnel inlet. 
Fig. I Different views of the inlet model. 
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a) View looking from outside into the tunnel inlet 



Fig, 2 Various perspective views of the inlet model of the NFAC Open Return Wind Tunnel 
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c) View looking from behind at the inlet model 


Fig 2 Concluded 
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a) Vector plot showing general behavior of the flow 


b) Enlarged view of velocity vectors near the inlet cowl 


c) Smoke flow and the tufts showing the recirculation bubble 
immediately downstream of the cowl. 


Fig. 5 Velocity vectors in the z-plane. 


















TUNNEL FLOOR SCAN 


PREDICTIONS USING PANEL 
O METHOD AT y - 5 0 ft A 



Fig 9 Tunnel floor centerline pressure and velocity (resultant) 
variation with the downstream distance 


SIDE WALL PRESSURE VARIATION 
O EXPERIMENT AT Z - 3 5 ft 



Fig 10 Pressure variation on the side wall centerline with 
downstream distance x 
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Fig 13 Companson of crossflow plane velocity vectors at x = -27 ft a) Crossflow plane grid, b) predicted velocity vectors, 

c) velocity vectors from experiment 
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